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Abstract 



A simple formula for the ratio of the number of odd- and even-parity states 
as a function of temperature is derived. This formula is used to calculate the 
ratio of level densities of opposite parities as a function of excitation energy. 
We test the formula with quantum Monte Carlo shell model calculations in the 
(p/ + (7g/2)-shell. The formula describes well the transition from low excitation 
energies where a single parity dominates to high excitations where the two 
densities are equal. 
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Parity is a fundamental property of nuclear levels, and its statistical distribution is 
important for describing parity-violating processes and neutron-capture reactions. Most 
theoretical models for level densities are based on the Fermi gas model Shell corrections 
and correlations due to residual interactions are included empirically. An empirical modifi- 
cation of the Fermi gas formula - the backshifted Bethe formula (BBF) - was successful in 
fitting many experimental level densities by adjusting both the single-particle level density 
parameter and the backshift parameter 0. Only limited data are available for the parity 
dependence of level densities since the neutron p-wave resonances are much weaker than the 
s-wave resonances at low energies and more difficult to measure. Ericson argued that 
the excitation of a relatively small number of single-particle levels with opposite parity can 
lead to an equal number of even- and odd-parity many-particle densities. The assumption of 
equal densities of opposite parities in the neutron resonance region is commonly accepted Q 
and used in the calculations of neutron-capture rates for s and r processes in nucleosynthesis 
0. Yet various theoretical studies [0 as well as analysis of experimental data ^ indicate 
that level densities can have a significant parity dependence. 

Parity properties can in principle be calculated within the interacting shell model, the 
basic theory of nuclear structure. However the calculation of level densities in the shell model 
requires large model spaces that are often beyond the reach of conventional diagonalization 
methods. Such methods are presently limited to A ^ 50 [pip!0| (in the p/-shell). Recently, 
quantum Monte Carlo methods |lTT| , p!2| were used to calculate total and parity-projected level 
densities in the framework of the interacting shell model. The methods were applied 



to nuclei in the iron-to-germanium region using the complete {pf + (79/2)-shell model space. 
The total level densities were found to be in good agreement with the experimental level 
densities, and significant parity dependence was found for A ^ 65. 

The Monte Carlo calculations of the parity-projected level densities accurately take into 
account shell effects and correlations due to the residual two-body interactions, but they are 
computationally intensive. In this paper we derive a simple formula for calculating the ratio 
of the number of odd- and even-parity states as a function of temperature. The formula 



is applied to nuclei in the iron region and compared with the Monte Carlo calculations. It 
reproduces well the crossover from low temperatures, where one parity dominates, to higher 
temperatures, where both densities become equal. Using the BBF for the total level density, 
the results of the model can be converted to a ratio of parity-projected level densities at 
fixed excitation energies. 

The Monte Carlo approach is based on the Hubbard- Stratonovich representation of the 
many-body imaginary-time propagator, e~^^ = / D[a]G{a)Ua, where G{a) is a Gaussian 
weight and 11^ is a one-body propagator that describes non-interacting nucleons moving 
in fluctuating time-dependent fields cr(r). The canonical thermal expectation value of an 
observable O can be written as {0)a = J -D[cr]G'(o")Tr^(0[/o-)/ / D[(r]G{(T)TTAUo-, where Tr^ 
denotes a trace in the subspace of A particles |14[. The integrand is easily calculated by 
matrix algebra in the single-particle space, and the multi-dimensional integral over the a 
fields is evaluated by the Monte Carlo methods. 

Parity-projected level densities were calculated in the Monte Carlo method using the 
projectors P± = (1 ± P)/2, where P is the parity operator [O. For even-even nuclei, the 



odd-parity level density is found to have large statistical Monte Carlo errors at lower energies 
(even for good-sign interactions) because of a sign problem introduced by the projection on 
odd-parity states. However, the sign problem does not affect the odd-even ratio of partition 
functions with the estimator 
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w. 
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Here ({a) = TiAUa and Cp{'^) = TrA(-Pf/cr)- In (|T]) we have used the notation {X„)w = 
J D[a]W{a)X^/ J D[a]W{a), where W{(t) = G{a)TTAU^. For a good-sign interaction 
and an even-even nucleus, W is positive definite. In the Monte Carlo method we sam- 
ple the fields a according to W{a) and then estimate {X^)\y ^ J^k^a-k/^^ where Ufc 
are M samples of the fields. At low temperatures {Cp{<^)/C{(^))w ~ 1, giving rise to a 
sign problem for the odd-parity states. The statistical Monte Carlo errors of Z±/Z are 
strongly correlated (since Z^/Z + Z^/Z = 1) and the error estimate of the ratio in (0) is 



A(Z_/Z+)^(Z+/^)-^A(Z+/Z). 

We have calculated the ratio Z^jZ^ for nuclei in the iron region using the complete 



(p/ + fl'9/2)-shell and the good-sign interaction of Ref. |jT3|. This interaction properly includes 



the dominant collective components of realistic effective nuclear interactions |T5|. In Fig. |1| 
we show the ratio of odd- to even-parity states as a function of inverse temperature (5 for 
three nuclei in the iron region: ^^Fe, ^°Ni and ^'^Zn. In all three cases we observe a transition 
from mostly even-parity states at low temperatures to an equal number of opposite parity 
states at high temperatures. However the crossover depends on the nucleus. For example, 
the crossover occurs at lower temperatures for ^^Zn than for ^^Fe. 

This observed parity dependence can be explained quantitatively by a simple model. We 
divide the single-particle levels into two groups of even- and odd-parities, and denote the 
group having the smaller average occupation by vr. If the particles occupy the single-particle 
states independently and randomly, we expect the distribution of occupancies n of the tt 
parity group to be Poisson: 

Pin) = L-e-f . (2) 

Here / is the average occupancy of orbitals with parity vr, which depends on temperature. 
We have assumed that the total number of the particles is sufficiently large compared with 
/• 

For a nucleus with even A, an odd-parity many-particle state is obtained when n is odd. 
The probability to have an odd-parity state is thus P_ = Yl'^^^ P{f^) = sinh/, while the 
probability for an even-parity state is = P{n) = e~-^cosh/. The ratio is then 

where we have identified P± = Z±/Z in terms of the partition functions Zj-{i3) of the 
even/odd parity states and the total partition function Z{(3). 

The argument leading to Eq. (^) is easily extended to the case where protons and 
neutrons are treated separately. We denote by Pp{np) and Pnirin) the respective Pois- 



son distributions for protons and neutrons with average occupancies fp and /„ of single- 
particle states with parity tt. A many-particle odd-parity state results if protons and neu- 
trons have overall parities (+, — ) or (—,+), respectively, i.e., P_ = Pp+P„_ + Pp_P„+ = 
e'^-^p^-^"-* (cosh /p sinh /„ + sinh /p cosh /„) (for an even-even nucleus). Similarly we find 
P+ = Pp+P„,_|_ + Pp_P„_ = e^*^-^p+-^"^ (cosh /p cosh + sinh /p sinh /„). The parity ratio is 
then given by Eq. (||) but with / = fp + fn- Furthermore, the convoluted distribution 
P{n) = Pp{np)Pn{nn) of finding n = Up + Un nucleons in orbitals with parity tt is by 

np+n„=n 

itself a Poisson distribution with f = fp + fn- 

Above the pairing transition temperature and in the independent particle model, / = (n) 
is evaluated from the Fermi-Dirac distribution {n) = J^aeni^ + exp[/?(ea — where the 

sum is over all orbitals of parity vr, and the chemical potential fi is determined from the 
total number of particles (in practice we use different chemical potentials for protons and for 
neutrons). To mimic interaction effects we use a single-particle spectrum that corresponds 
to an axially deformed Woods-Saxon potential. For axial deformations, parity is a good 
quantum number, and we can estimate (n) by summing over the Fermi-Dirac occupations 
of all deformed orbitals with parity vr. 

To test how well the Poisson distribution (^ describes the distribution P{n) of the 
occupation of the single-particle states with parity tt, we compare with the Monte Carlo 
results. The probability to find n particles in single-particle states with parity tt is 

_ Tr^ [e-'^'Sih - n)] _ / Tr ^ [UJ ~ n)] \ 
^ Tr^e-^^ -\ Tr^f/. " 

The quantity inside the brackets of Eq. (HI) is calculated from a double projection on particle 



number A and occupation number n of states with parity tt 

1 N 

Tta [UJin - n)] = Yl E e'^'^-'^e-*'^^-" det (l + e"^-e"^'^'-V„) , (5) 

where /,r is a diagonal matrix with diagonal elements 1 for each orbital of parity vr and 
otherwise. The quadrature points are given by (pm = 2TTm/N and = 2vrA;/(A^7r + 1), where 
is the total number of single-particle states and N^^ is the number of orbitals with parity 
vr. Uo- is the N x N matrix representing the propagator Ua in the single-particle space. 



For nuclei in the iron region the occupation of the even-parity orbital (^9/2 is relatively 
small, and we choose n to denote the occupation of the (79/2 states. Using the complete 
ipf + 5'9/2)-shell we calculated the distributions P{n) from Eqs. and (^. The results 
are shown in Fig. ^for ^^Fe, ^''Ni and ^^Zn at several temperatures (solid circles). The solid 
lines are the Poisson distributions (0) with / taken to be the average occupation of the (79/2 
orbital (calculated in the Monte Carlo). At high temperatures the microscopic distributions 
are well described by the Poisson distribution. However, for temperatures T ^ 1 MeV we 
observe deviations describing the enhancement of P{n) for even n and the suppression for 
odd n due to pairing effects. 

At lower temperatures T < 1 MeV, it is necessary to take into account pairing effects. 
We can still use Eq. (^, but now with quasi-particles. Consequently, Eq. (Q) is applicable 
where n is replaced by the number of quasi-particles with parity vr, and / in Eq. (^ is the 
average occupation of quasi-particle states with parity vr. 



As in the original BCS treatment |T6[, the occupation probabilities have a compo- 
nent from condensed pairs and a component fa from quasi-particles. Minimizing 
the free energy gives the quasi-particle occupation factor fa = + 6xp(/?-E'a)], where 



Ea = y {ea — + and the gap A and the chemical potential A (at finite T) satisfy 
self-consistency conditions. Since the condensed pair occupations play no role in the parity 
of the states, one should only use the quasi-particle / in Eq. (H) 

■^^?/"^£l + ex!)(/3E.)- 

We have applied the above model to determine the parity-ratio of levels for the three 
nuclei shown in Fig. [l|. The deformation parameter 6 (used to calculate the single-particle 
spectrum e^) is extracted from the experimental B{E2) values [jl^ for the 2+ 0+ transition 
using B{E2) = [{3/4:n)ZerlA^/^6Y/5. We find (using tq = 1.27 fm) 6 = 0.22,0.18 and 0.19 
for ^^Fe, ^°Ni and ^^Zn, respectively. The pairing gap A (at T = 0) is extracted from odd- 
even mass differences and is used to determine the pairing strength G by a BCS calculation. 
The total occupation / of the quasi-particle even-parity states (^ is shown by the solid lines 



on the left column of Fig. ^ Above the critical BCS temperature / coincides with (n). 
Below the BCS temperature the average number of particles with parity tt is different from 
/ and includes a contribution from condensed pairs 



H = E[/« + ^a(l-2/a)] = ^E 



2 



1 !fL-^. ^.(^^- 
1 — tanh 



(7) 



The calculated (n) is shown by the dashed lines in Fig. ^ The parity-ratio is 
calculated from Eqs. @ and (^) and shown by the solid lines on the right column of Fig. 
^. The model describes well the Monte Carlo results (symbols). 

Eq. d^) expresses the ratio of the number of odd- and even-parity levels at constant 
temperature. To calculate the odd-to-even level density ratio at constant excitation energy 
we use in addition the Lang and LeCouteur version |TB[ of the BBF for the total level density 
p{Ex). We calculate the total partition function from Z{f3) = e'^^^'' J p{Ex)e~^^'^ dE^-, where 
Egs is the ground state energy. Using Z^ + Z^ = Z and Eq. (§) for we can determine 

Z±{I3) = Zj (l+tanh^^ /) and calculate the thermal energies for even- and odd-parity states 
from E± = -ainZ±/5/3 = E + [rfln(l + tanh^^ f)/df]{df/dp). We can then calculate the 
canonical entropies and heat capacities from standard thermodynamic relations and find the 
parity-projected level densities p±{Ex). Fig. ^ shows the calculated ratio p_{Ex) / p+{Ex) in 
our model versus excitation energy (solid lines) for ^^Fe, ^°Ni and ^^Zn. The results compare 
well with the Monte Carlo calculations shown by error bars. 

In conclusion, we have derived a simple formula for the parity dependence of level densi- 
ties. The formula describes the crossover from low excitation energies where a single parity 
dominates to higher excitations where odd- and even-parity states have equal densities, and 
agrees well with shell model Monte Carlo calculations. 
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FIG. 1. Ratio of odd- to even-parity partition functions Z_ (/?) (/?) versus inverse tem- 
perature P for ^^Fe (circles), ^'^Ni (squares) and ^^Zn (diamonds), calculated in the Monte Carlo 
method according to Eq. ([l|). 
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FIG. 2. The distribution P{n) of the number of particles n occupying the even-parity 59/2 
orbitals. The sohd circles are the Monte Carlo results using (^). The solid lines are Poisson 
distributions (||) with the same values of (n) calculated in the Monte Carlo method. 
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FIG. 3. Left column: The solid lines are the occupation / of the quasi-particle even-parity levels 
@ versus P for ^^Fe, ^'^Ni and ^^Zn. Above the BCS temperature / coincides with the average 
occupation (n) of the even-parity states. The dashed lines show (n) below the BCS temperature 
(Eq. (0)) while the dotted-dashed lines are calculated from the Fermi-Dirac occupations. Right 
column: the ratio Z-jZj^ versus /3 calculated from using the occupations / shown on the left 
(solid, dashed and dotted-dashed lines). For comparison we show by symbols the Monte Carlo 
results. 
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Ex 

FIG. 4. The parity ratio p^{Ex) / p+{Ex) versus excitation energy Ex- The solid hnes are 
calculated from Eqs. (^) and (|6|) (see text), and the symbols with error bars are obtained in the 
Monte Carlo method of Ref. (shown with error bars). 
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